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The small r variation of the probability density P(r) for end-to-end separations of a -CH 2 CH 3 capped (- 
OCH2CH2-)„ oligomer in water is computed to be closely similar to the CH4 • • ■ CH4 potential of mean force 
under the same circumstances. Since the aqueous solution CH4 ■ ■ ■ CH4 potential of mean force is the natural 
physical definition of a primitive hydrophobic bond, the present result identifies an experimentally accessible 
circumstance for direct observation of a hydrophobic bond which has not been observed previously because of 
the low solubility of CH4 in water. The physical picture is that the soluble chain molecule carries the capping 
groups into aqueous solution, and permits them to find one another with reasonable frequency. Comparison 
with the corresponding results without the solvent shows that hydration of the solute oxygen atoms swells 
the chain molecule globule. This supports the view that the chain molecule globule might have a secondary 
effect on the hydrophobic interaction which is of first interest here. The volume of the chain molecule globule 
is important for comparing the probabilities with and without solvent because it characterizes the local 
concentration of capping groups. Study of other capping groups to enable X-ray and neutron diffraction 
measurements of P(r) is discussed. 



Hydrophobic interactions are central to super- 
molecular self-assembly in aqueous solutionf^^, includ- 
ing the folding of soluble globular proteins™ A sur- 
prizing development of recent years is that for entropy- 
dominated hydrophobic solubilities we now have statis- 
tical mechanical theories that are fully defensible on a 
molecular scale, exploiting all available molecular-scale 
dataPHUl In contrast, our understanding of the statis- 
tical solvent-induced forces between neighboring small 
hydrocarbon molecules in water — hydrophobic interac- 
tions — has not experienced t he co rrespondingly conclu- 
sive progress^ despite intense^ESl computational effort. 

The principal impediment to progress in understand- 
ing hydrophobic interactions is the lack of an accessible 
direct observation of a primitive hydrophobic bond, such 
as the CH4 ■ ■ ■ CH4 potential of mean force. Because 
simple hydrophobic species, such as CH4, are only spar- 
ingly soluble in water, direct observations of primitive 
hydrophobic interactions are difficult to achieve; instead 
the influence of hydrophobic interactions is inferred in 
more complex systems — protein folding being the fore- 
most example!^ Where limited thermodynamic evalu- 
ation of hydrophobic interactions has been experimen- 
tally achieved for soluble, but more co mplex hydrophobic 
species with higher aqueous solubilities j 51 l 52 l those exper- 
iments set-off a substantial modeling effort that has not 
untangled this complicated issuePS 

Here we show (FIG. [I]) that the probability density 
P(r) of end-to-end separations of a -CH 2 CH 3 capped (- 
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OCH2CH2-)„ oligomer in water exhibits a distinct hy- 
drophobic bond between the two end C-atoms. These 
atom-pair correlations are intrinsically measurable by 
X-ray and neutron diffraction.^ A natural view of the 
present case is that the chain molecule carries this hy- 
drophobic pair into solution and permits them to find 
one another with reasonable frequency. The underlying 
assumption is that the hydrophobic interaction is suf- 
ficiently local that the effects of the supporting chain 
molecule are secondary. This is a reasonable assump- 
tion that can be experimentally tested by variation of 
the chain length and capping groups. Nonetheless, the 
present realization is consistent with the view that hy- 
drophobic interactions are typically expressed in the con- 
text of other effects in micelles, membranes, and the 
structure of soluble proteins. 

Two distinct sets of calculations were combined to 
obtain the present results: firstly, parallel tempering^ 
to establish the overall structure of P(r) and, secondly, 
windowing^ to achieve satisfactory spatial resolution in 
the interesting loop-closure regime, r < 0.7 nm. Direct 
observation of the contact feature would rest on a few 
percent of the parallel-tempering data set, i.e., about 
2% of the molecules are in loop-closure configurations. 
The overall parallel-tempering results permit evaluation 
of the undetermined scale factor (zero of the potential of 
the average forces) for the windowing results. 

The shape of P(r) in the loop-closure region is strik- 
ingly similar to pr edicted pair distributions for model 
inert gases in water PH 61 I 62 I On the basis of the stratified 
evaluation of P(r), the most probable -CH3 • • • CH3- 
hydrophobic bond length is about 0.4 nm. The slight 
shoulder on the large-r side of the principal maximum 
of P(r) is a remnant of a RISM cusp^ and reflects the 
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FIG. 1. Probability density P(r) for radial displacement of 
end-C atoms of CHaCHaOtCHaCHaOjaoCIlaCHs in aque- 
ous solution, with the cumulative distribution J Q r P(x)Aitx 2 &x 
(dashed line) . The dots were obtained by parallel tempering^ 
over times of 20 ns/replica with 32 temperatures spanning 
256-550 K, using AMBER 1(P with the same molecule num- 
bers and cubical volume V = 32.254 nm 3 for each replica. 
Monte Carlo swaps of configurations had an overall success 
rate of 23%. The SPC/E model treated the 1000 water 
moleculeip^, and an extended atom model with the Gener- 
alized Amber Force Field (GAFF) was used for the single 
chain molecule.^ The solid line is a stratified recalculation 
of the loop-closure feature, adopting a harmonic potential for 
the radial displacement coordinate covering the range of 0.3- 
1.0 nm uniformly with 15 windows, then reconstituting P(r) 
on that range with a weighted histogram methodP^^ This 
second set of calculations lasted 20 ns/window. The undeter- 
mined multiplicative constant in the stratified calculation of 
P(r) is adjusted to match the direct observation from the par- 
allel tempering in the region of the first minimum. The left 
vertical axis is non-dimensionalized with the observed (r 2 ) 1//2 
= 1.56 nm. 



ethyl end-capping. 

A direct comparison with results for CH4 • • • CH4 ob- 
tained with the same models and stratification methods 
(FIG. [2]) shows good agreement in the placement of the 
principal peak, and the large-r shoulder of the principal 
peak is more obvious. Contact pairing configurations are 
more prominent than solvent-separated ones in both of 
these cases, supporting a standard view of this hydropho- 
bic interaction. A separate evaluation (not shown) of the 
solubility of CH^aq) with these models compares well 
with experiment. This substantiates the view that this 
model solute naturally represents hydrophobic methane. 
These probability densities differ in outer shells due to 
detailed differences between the solvent and the chain 
molecule medium. 

It is natural also to compare these pdfs with what 
would be obtained if the solvent were absent (FIG. [3 1. 
The comparison shows that water swells this soluble 
chain molecule (FIG. [IJ), as is the case also with aqueous 



FIG. 2. Comparison of the probability density of FIG.[TJ(solid 
curve) with results from the CH4 • • • CH4 potential of mean 
force obtained with the same models and stratification meth- 
ods (dashed curve) as FIG. [TJ The latter result is in good 
agreement with a previous simulation evaluation for a similar 
model.^ The multiplicative constant required for this com- 
parison was adjusted informally for agreement in the region 
of the maxima of these two results. 



dilution of a PEO meltP^ This supports the view that 
the chain molecule globule might have a secondary effect 
on the hydrophobic interaction which is of first interest 
here. 

Note that the boundary between the low-extension 
r < (r 2 ) region and a high-extension region is def- 
inite (FIG. [l]). On the basis of the similarity with the re- 
sults for FIG. [2j the end-cap pair appears to be shielded 
within a uniform fluid environment in the low-extension 
region. Furthermore, end-caps are more concentrated in 
the smaller- volume polymer globule of FIG. [3] than in 
the larger- volume polymer globule of FIG. [TJ The non- 
dimensionalized probability densities incorporate that 
distinction without which the maximum probabilities 
would differ by nearly a factor of three. This point is 
consistent with the observation from FIG. [2] that if the 
probability densities there match roughly at large-r, then 
the maximum values also approximately match. 

The high-extension r > (r 2 ) 1 / 2 tail of the probabil- 
ity density (FIG. [TJ is reasonably described by a linear- 
response behavior 



dr 



In P(r) 



br 



(1) 



A Gaussian model for P(r) at high-extension would not 
be centered on the origin, however, and in that sense a 
traditional Gaussian model would be unsatisfactory in 
the high-extension region. 

Simple characterizations of hydrophobic effects, i.e. 
whether they are attractive or repulsive, depend on the 
specific properties examined and the comparisons made. 
For example, well-developed theoretical analyses-^ show 
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FIG. 3. As in FIG. [T] but without water present. The left 
vertical axis is non-dimensionalized with the observed (r 2 ) 1//2 
= 1.06 nm. 



that hydrophobic interactions can be repulsive for the os- 
motic second virial coefficient. Other comparisons illumi- 
nate different aspects of hydrophobic interactions. Some- 
times hydrophobic interactions are judged by adopting 
another solvent that provides a natural comparison to 
the case of the water medium!^ Sometimes it is most 
direct to judge hydrophobic interactions by comparison 
with standardly hydrophilic solutes, polar, H-bonding, or 
ionic species in water. We propose that identification of 
an experimentally accessible case that also permits de- 
tailed molecular theory and computation should assist in 
resolving such alternatives. 

Though this hydrophobic interaction is intrinsically 
measurable, analysis of X-ray and neutron diffraction ex- 
periments on such systems will require specific and thor- 
ough support from molecular simulations, as is the cur- 
rent practice! 53 * 66 ^ Isotopic substitution and labeling will 
be essent ial. Consideration of fluorinated caps such as - 
CF j67 | 68 | ghould make this hydrophobic interaction more 
prominent yet.^ 

The chain molecules considered here have broad tech- 
nological interest because of their biocompatibilities. 69 
They are also intrinsic to the dispersant materials used in 
response to oil spills. 70 Small angle neutron scattering^ 
and fluorescence and light-scattering studied have 
shown that the aqueous solution interactions of (- 
OCH2CH2-)„ polymers are sensitive to the end-capping 
of the chains. The importance of end-effects is also sup- 
ported by the sensitivity o f solu tion phase diagrams to 
the (-OCH 2 CH 2 -) n lengths J^ 17 -^ These observations sug- 
gest hydrophobic bonding of the -CH3 caps, and that ma- 
nipulation of capping groups might help in understand- 
ing hydrophobic interactions on a molecular scale^ by 
exhibiting a localized hydrophobic bond. 

We thank D. Asthagiri and M. Berkowitz for helpful 
discussions. 
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